Spatial decoherence near metallic surfaces 
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We present a first-principles derivation of spatial atomic-sublevel decoherence near dielectric and 
metallic surfaces. The theory is based on the electromagnetic-field quantization in absorbing dielec- 
tric media. We derive an expression for the time- variation of the off-diagonal matrix element of the 
atomic density matrix for arbitrarily shaped substrates. For planar multilayered substrates we find 
that for small lateral separations of the atom's possible positions the spatial coherence decreases 
quadratically with the separation and inversely to the squared atom-surface distance. 
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I. INTRODUCTION 

New physical models for quantum information process- 
ing and quantum computation have been inspired re- 
cently by the experimental achievements in trapping and 
controlling ultracold neutral atoms 0, 0, 0, US ■ The 
first experimental step to achieve a physical realization 
of a quantum computer with neutral atoms is to confine 
them on a definite region in the space. The creation of 
microscopic guides and traps for neutral atoms moving 
close to surfaces is possible using nanofabricated struc- 
tures that either carry currents or are based on perma- 
nent magnetized films. The idea at the base of atom 
chips has been put forward by Frisch and Segre [|| who 
realized that, when a homogeneous magnetic field ('bias 
field') is superimposed with the field created by a cur- 
rent flowing through a wire, the magnetic field vanishes 
on a line parallel to the current which can trap atoms in 
low- field seeking magnetic hyperfine sublevels. 

One of the main requirements for a qubit is to be well 
isolated from a noisy environment to avoid decoherence, 
namely the destruction of quantum superpositions due 
to the coupling of the atom cloud to the noisy chip en- 
vironment. Although neutral atoms are considered good 
candidates as quantum systems since they have a small 
coupling to the environment, they still suffer from loss 
and decoherence. When atoms are trapped in atom chips, 
they are held close to the material surfaces. The small 
separation between the cold atom cloud and the macro- 
scopic environment (usually at room temperature) raises 
the question of how strong the energy exchange will be, 
and which limit of atom confinement and height above 
the surface can ultimately be reached. Thermal fluctua- 
tions induce noise currents in the materials the trap 
is made of, and fluctuations of the electromagnetic field 
are produced in the conducting body. Such fluctuating 
fields can be strong enough for an atom close to the sur- 
face to drive rf magnetic dipole transitions that flip its 
spin causing either its loss or decoherence of its quantum 
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state. 

In [E HO, dH O HE Ull , atom loss due to thermally 
driven spin flips has been widely investigated and sev- 
eral experiments have confirmed the theoretical findings 
(la. fl6t Il7| . In this article we examine the influence of 
thermally-induced spin flips on the coherence properties 
of atomic spatial superposition states. Such coherent su- 
perpositions can be thought of being created by tunneling 
through a shallow potential barrier in either a double- well 
potential or, more generally, an optical lattice structure 
[l8| . The study of the latter has been received much at- 
tention over the recent years for its potential application 
in quantum information processing (see, e.g. |19l |20(). 
The derivation is carried out within the framework of 
the quantum electrodynamic theory for electromagnetic 
fields in dielectric media H3, E3, E3, H3, Hg, Hg, [23, H^l 
which yields a first-principle description of the decoher- 
ence properties of spatial atomic superposition states. 

This work is organized as follows: Sec. II introduces 
the basic notions of a quantized electromagnetic field in a 
dielectric medium. In Sec. Ill the density matrix of the 
atom is obtained in the presence of a fluctuating mag- 
netic field and an expression for the spatial coherence is 
derived. We focus on a particular substrate geometry, a 
planar ly multilayered structured, in Sec. IIVI for which 
the dyadic Green function is explicitly known. 



II. BASIC EQUATIONS 

It is well known that the quantum statistical proper- 
ties of electromagnetic fields and their interactions with 
atomic systems can be strongly influenced by the pres- 
ence of dielectric bodies. In the present context it is 
useful to formulate quantum electrodynamics 
a dielectric-matter background 

I27I |2S[ . The interaction between atomic systems and 
the electromagnetic field is typically treated in terms 
of the polarization and magnetization associated with 
the atomic charges. Let us restrict our attention to an 
isotropic but arbitrarily inhomogeneous medium whose 
polarization responds linearly and locally to the electric 
field. Causality and the dissipation-fluctuation theorem 
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[29j then require that 

oo 

P(r,t)=e / dT X {r,t)E(r,t-T)+P N (r,t), (1) 
o 

where x(r, t) is the dielectric susceptibility (in the time 
domain) and Pjv(r, t) is the noise polarization associated 
with dissipative processes in the dielectric medium. 

Using Maxwell's equations in Fourier space, we find 
that E(r,w) obeys the Helmholtz equation 

uj 2 

VxVxE(r,w) 5-e(r,w)E(r,w) = LJ 2 fj, P N (r, u), (2) 

tr 

where the complex permittivity, e(r, ui) — er{v,uj) + 
iei(r,uj), is defined by 



e(r,w) = 1+ / dre Mt x(r,r). 



(3) 



The solution to Eq. can then be written as 

E(r, uj) = u; 2 /io / d 3 r'G(r, r', ^P^r', «), (4) 

where the Green tensor G(r, r', w) is a second rank tensor 
that has to be determined from the partial differential 
equation 

uj 2 

V x V x G(r, r',ui)- -=-e(r, w)G(r, r', w) = <5(r - r')U, 
&■ 

(5) 

where t/ is the unit dyad. An important consequence of 
the differential equation JSJ) is the integral relation [2{| 

f uj 2 

/ d 3 s— si(s,uj)G(r,s,uj)G + (r',s,uj) = ImG(r,r',u). 

Quantization of this theory then proceeds in the usual 
way [2^] . First, a factor is split off from the (classical) 
noise polarization, 



Pjv(r, uj) = i\j — £/(r, uj) f (r, uj). 



(7) 



One then identifies the dynamical variables f (r, uj) as the 
fundamental 5 correlated Gaussian random process and, 
upon quantization, replaces them by the operator-valued 
bosonic vector field f (r, uj) satisfying the equal-time com- 
mutation relations f (r, uj), f'(r', uj') — S(r — r')5(uj — 

uj')U . The Hamiltonian of the system composed of elec- 
tromagnetic field and absorbing matter is 



00 

H F = J d 3 r J dujhuj&(r,uj)i(r,uj). 



(8) 



The electromagnetic field operators can now be ob- 
tained in the Schrodinger picture as 



E(r) = / rfcjE(r,w) + H.c, (9) 


E( r , w )=i A /JL^ f £ fVV£i(r',w)G(r ) r»f(r» 
V tteq c z J 

(10) 

and, using Faraday's law, 

B(r,w) = (ia;) -1 V X E(r,w). (11) 

An important feature of this theory is that it repro- 
duces the correct form of the fluctuation-dissipation the- 
orem. Let the system of electromagnetic field and ab- 
sorbing matter be in thermal equilibrium at some tem- 
perature T . Then the thermal correlation function of the 
dynamical variables at temperature T reads 

(f (r, w)f V, uj')) = (n th + 1)8(t - v')S(uj - uj')U, (12) 

with the mean thermal photon number at frequency ui 

1 



n t h 



e huj/k B T _ -y ' 



(13) 



From Eqs. (JTUJl — (H2J), together with Eq. ©, it follows 
that the thermal expectation value of an anti-normally 
ordered product of magnetic field operators can be writ- 
ten as 



(B(r, W )Bt(r>')> = 

%0 n 



-Im 



V x G(r,r',uj) x V (n th + l)S(u - uj'). (14) 



Such a quantization model provides a valid description 
of electromagnetic field in absorb ing dielectric materials. 
In fact, it has been shown in 0iEl| that the equal-time 
basic commutation relations of QED are preserved. The 
electromagnetic field is expressed in terms of the classi- 
cal Green tensor satisfying the Helmholtz equation J5J, 
and the continuum of the bosonic field variables i(r,uj). 
All the information about the dielectric matter is con- 
tained in the Green tensor via the permittivity e(r,w). 
For metals at low frequencies, the permittivity can be 
approximated by the well-known Drude relation 



e(w) 



2ic 2 



2 S 2 



(15) 



with the skin depth 5. Although such a relation is not 
strictly consistent with causality as it has recently been 
pointed out [30j, it can be assumed to be valid in a re- 
stricted frequency interval. 

At this point it is necessary to point out the limita- 
tions of the quantization scheme presented above. Note 
that the form of the polarization, Eq. QJ, is valid only for 
strictly locally responding materials. That is to say, we 
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assume that the elementary dipoles that give rise to the 
polarization are essentially fixed in space. Certainly, for 
metals which can alternatively described by a conductiv- 
ity, this is not true as charge carriers can move around 
freely for considerable distances. However, the locality 
assumption can be upheld in situations in which the mean 
free path length is much shorter than all the other length 
scales in the system under consideration. While this is 
certainly true for ordinary metals at room temperature 
and geometric length scales of several micrometers, we 
do expect corrections due to spatially nonlocal response 
(the anomalous skin effect) for metals or superconductors 
at very low temperatures as considered in 0, ^| . 



where the operator of the magnetic induction is given by 
Eq. (JTTJ, together with Eqs. |J3J and iJTUjl. The magnetic 
moment operator in Eq. I|17[) associated with a transition 
\i) — > |/) can be written as fi = fi\i)(f\ + h.c. Since we 
assume the atom to be cooled into its electronic ground 
state, there is no contribution of the angular momentum. 
Furthermore, since the nuclear magnetic moment can be 
neglected because of the ratio of the electron mass to 
the mass of the nucleus (see the discussions in [9I IT^V 
the magnetic moment vector is just proportional to the 
expectation value of the electronic spin operator, 

» = gstiB(i\S\f), (18) 



III. SPATIAL DECOHERENCE 

Let us suppose we had an atom in one of two adjacent 
sites of an optical lattice. The tunneling interaction al- 
lows the atom's wave function to coherently spread over 
the neighboring site [3l| where its state can be written 
in the occupation-number basis as 



m = o)) A 



V2 



(|1,0) + |0,1)). 



(16) 



We take the time at which the equal superposition has 
been established to be t — and assume for simplic- 
ity that no tunneling occurs at later times, at least not 
at timescales shorter than the decoherence time. This 
means that we imagine the tunneling interaction being 
frozen over a certain time period. This assumption is jus- 
tified when considering proposals in which spatial atomic 
locations are used to encode quantum information. 

Atoms that are held close to microstructured surfaces 
experience fluctuations of the electromagnetic field due 
to absorption in the substrate material. In the case of 
a magnetic trap the atom is subject to a constant mag- 
netic field with strength Bq in the center of the trap. 
The magnetic sublevels are split due to the Zeeman ef- 
fect by the Larmor frequency lol — gsfJ-BBo/h. A subset 
of these magnetic sublevels feel an attractive potential 
towards regions of low magnetic field. In the experiment 
reported in ^{| 87 Rb atoms are initially pumped into 
the hyperfine state \F, mp) = 1 2, 2) in which they are 
trapped. However, due to absorption in the surface ma- 
terial and the resulting quantum fluctuations, fluctuat- 
ing magnetic fields cause the atoms to evolve into states 
with lower magnetic quantum number mp. In sufficiently 
tight magnetic traps, also atoms in the \F,itif) — |2, 1) 
state are trapped. Spin flips to even lower magnetic sub- 
levels cause the atoms to be expelled from the trap. In 
this case, spatial decoherence is no more a matter of in- 
terest. Hence, it is sufficient to treat the atomic system 
in a two-level approximation. 

We focus on the Zeeman coupling of the atomic mag- 
netic moment to a fluctuating field represented by the 
Hamiltonian 



Hp. 



-A-B(r A ) 



(17) 



where fis denotes the Bohr magneton, and gs ~ 2 the 
electron's g-factor. Inserting Eq. (|llfl into Eq. (|17fl . the 
Zeeman Hamiltonian can be written in the rotating-wave 
approximation as |T^ | 



H 7 = - 



^Bgs 



-Hb9s 



{f\S q \i)£}B q (v A )+h.c. 



(f\S q \i) J daj^^J -^-e qpj d p / <l A s 





X yj £/(s, Lj)Gji(r A ,S, U))fi(s, + h.i 



(19) 



where 5 = \ f)(i\ denotes the atomic spin lowering opera- 
tor. Finally, the free atomic Hamiltonian can be written 
in the two-level approximation used above as 



H A = ^aL = \tkJ A (m\ - l/X/l), (20) 



where the £ obey the commutation rules , £ z ] = =f£^ ■ 
In order to analyze how this magnetic noise influences 
the coherence of the state of our atom, we rewrite the 
initial atomic state as 



1^ 



1 

V2 



l*2». 



(21) 



where the labels 1,2 refer to the occupied site. Let us 
consider a system composed of the two-level atom and 
a fluctuating magnetic field initially in the vacuum state 
|0), so that the total state of the atom-field system reads 



1 



V2 



(|il,0) + |i 2 ,0». 



(22) 



The Hamiltonian describing the evolution of the com- 
bined system is given by the sum of the three Hamilto- 
nians H — Hp + H& + H z , where Hp, and H Zl and Ha 
are given by Eqs. ©, (|2U|) . and (jT§|) . respectively. The 
system wave function at a certain time t can be written 
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as m 



\iiAF(t))=C ll (t)e- t ^ 2 \t 1 ,0)+C l2 (t)e-^ 2 \t 2 ,0) 

oo 



+ J d 3 r I ^C /2 , m (r, W ,t)e-' i( "-^ /2) *|/2,l m (r,w)), 
o 

(23) 

where |0) and |l m (r, ui)} denote the electromagnetic field 
vacuum and single-excitation states, respectively. The 
Schrodinger equation ihd t \^AF{t)) = H\ip AF (t)) yields 
(a = 1,2) 



C ia (t) = J^L{f\S q \i) I ,/ ! r / ,U 



xcoe 



—i{uj—uJA)t 



\/ei(r,u))e qp jdpG jm (r a , r, u>) 



xCf atfn (r,u,t), 



(24) 



C fa , m (T,u,t) = -^B={i\S q \f)cje^-^ 

x y/ei(r, uj)e qp jd p G* jm (r a , r, uj)C ia (t). (25) 

We now substitute the result of formal integration of 
C/ a . m (r,w,t) with the condition C/ a . m (r, oj, 0) = into 
Ci a (t) , make use of the integral relation @) , and obtain 

C\ a (t) = f dt' K a (t - t')C ia (t'), (26) 
Jo 

where the integral kernel is 

K a (t-t') = -^^{f\S q \i) (i\S k \f) 

oo 

x J 4, e -'(^^)('-0 Im 
o 



V x G(r a ,r a ,uj) x V 



(/ k 

(27) 



We integrate both sides of Eq. (|26|l over t, and change 
the order of integrations on the right-hand side we derive 



and the initial condition Ci„(0) = 1. When the Markov 
approximation applies, i.e., when in coarse grained de- 
scription of the atomic motion memory effects are disre- 
garded, we may let 

'- -> -ttS(uj - oj a ) + iV 

t(u> - uj a ) 

Defining the coefficients 



ui — ui A 

(30) 



r, = 2(^f^] (f\s q \i){i\s k \f) 



xlm 



V x G{v a ,v a ,u A ) x V 



qk 



(31) 



and 



Im 



xV J duj- 

o 



V x G(r a ,r a ,uj A ) x V 



<jk 



OJ — uja 



(32) 



we can write K a (t — t') = — ^T a + i5uj a . We finally obtain 
for the time evolution of the coefficients Ci a (t) 



C ia (t) = exp 



-T a + iSu} a ) t 



(33) 



The coefficients T a and 8uj a defined in Eqs. (|31|l and l)32|) 
represent the spin flip rate and the line shift, respec- 
tively, and have been derived in a similar fashion in |l2j . 
The spin flip lifetimes 1 /T„ have already been subject of 
major theoretical H H Efl IH El El H and experi- 
mental 0, 0, investigations which will not repeated 
here. In what follows, we will assume that the line shift 
Suj a caused by the interaction with the quantized elec- 
tromagnetic field is negligible. This can be seen as fol- 
lows. The Green function appearing in Eq. as well 
as the Fourier transform of the permittivity in Eq 
plays the role of a response function and so it satisfies 
the Kramers-Kronig relations for a complex- valued func- 
tion g(uj) = Re[g(w)] + ilm[g(u>)} 0, 



Re[g(u) 



duj 



,Im[ 5 (w)] 



(34) 



C ia (*) - C ia (0) = J dt' K a {t - t')C ia (£') (28) 
o 

with 

Ka(t-t') = ^^-(f\S q \i) (i\S k \f) 



Z 2 TT£f)h 



e -i(o)-o)A)(t-t') _ I 

x / duj - ; Im 



i(u - lu a ) 



V x G(r a ,r Q ,w) x V 



(29) 



Im[ff(6 



= V 

TV 



, Re [g(u 



ui — UJ 



(35) 



The lower limit of the integral in Eq. (|32|l can be ex- 
tended to — oo with little error as the integrand is peaked 
around u A . Hence, Eq. Il32[l can be rewritten as 



5w a 



C 2 £oH 



(f\S q \i) (i\S k \f) 



xRc 



V x G{v ai r ai Lu A ) x V 



- qk 



(36) 
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As we will see later, the line shift is of the same order of 
magnitude as t he sp in flip rate. For typ ical exp erimental 
realizations, 0, EDt [HI IS E3, E3, EH EH , this will 
be in the sub-Hz range. This means that 5tv a can be 
neglected as it is extremely small when compared to the 
spin flip transition frequency. 

Now substituting Eq. l-i-'il) into the expression for 
Cf a>m (r,uj,t), Eq. H25|) . we find the formal solution 



C/ , m (r,w,*) 



[i\S q \f)wy/ £/(r, aj)e qpj d p 



o 



(37) 



In order to find how the off-diagonal elements of the den- 
sity matrix decay, we trace the atomic density matrix 
over the field and obtain 

g A (t) = (0\g AF (t)\0) 

oo 

+ Y^J d3r J diu(U(r,u})\QAF\U(r,uj)) 



1 ( Pu(t) pi 2 (t) 

2 I Pl2(t) P22(t) 



(38) 



where the matrix elements Qij of the density matrix have 
to be calculated from 



Qn(t) 

022 (*) 



OO 

\C h (t)\ 2 +Yl J d3r J du\C fum (r,uj,t)\ 2 , 



(39) 

OO 

\C l2 (t)\ 2 +Y, J d3r J du>\C f2!m (r,oj,t)\ 2 , 
1 o 

(40) 

OO 

, m (r,w,£) Cj 2 m (r, u,t), 

o 

(41) 



First, it can be checked that the diagonal elements Qn(t) 
and Q22(t) are properly normalized to Qu{t) = Q 2 2(t) = 1 
by inserting Eqs. I|33[) and l|37|l together with Eq. JHJ into 
Eqs. I|39|l and l]4Up. respectively. Thus, as a consistency 
check we find that Tr[gyi] = 1. We can then calculate the 
off-diagonal elements of the density matrix as 



012 (*) = ' 



+ 2(1- 



' c 2 eoh 



Im 



x{i\S q \f)(f\S k \i) 



V x G(r 2 ,ri, uj a ) x V 



kq 



where ]?i2 = (ri + T2)/2 is the arithmetic mean of the 
spin flip rates, Eq. H31|) . at both sites. Note that the 
Hermiticity of the density matrix QA(t) follows from the 
reciprocity theorem applied to the dyadic Green function 
which yields G(n, r 2 , uj a ) = G T (r 2 , ri, uj a ). 

Equation (|42|l constitutes the main result of our pa- 
per. It provides, via the Green function G(r 2 ,ri,uj A ), 
an elegant way to assess the loss of spatial coherence for 
arbitrarily shaped substrates. Recalling the expression 
for the fluctuation-dissipation theorem, Eq. i|14|) . it fol- 
lows that Eq. I|42J) can be rewritten as 



M<) = e- ri2t +(l-e- ri2t ) 
y WS q \f}(f\S k \z} MMr2^A)Bl(ru^)) 
(i\S q \f){f\S k \i) dw(4(ri,WA)B2(ri,w)) 



S {r 1 ,r 2 ,uj A ) 



(43) 



in terms of the magnetic cross-correlation tensor 
(B(r, £j)Bt(r', ui')). This means that the imaginary part 
of the (magnetic) Green function is proportional to the 
spatial coherence function of the fluctuating magnetic 
field ISl|3El3i. 



Note that, although the calculations have been per- 
formed for surfaces held at zero temperature, the exten- 
sion to finite temperatures is trivial. Indeed, it is seen 
from Eq. (|14fl that the spatial coherence functions as well 
as the spin-flip rates simply have to be multiplied by the 
factor (nth + 1) to account for thermal fluctuations. 

Equation (|42|l . or equivalently, Eq. consists of two 
parts. The first is a (spatially local) exponential decay 
that describes the effect of the transition from the initial 
spin state \i) to the final spin state |/). The second term 
is a (spatially nonlocal) non-exponential term which is 
proportional to the spatial coherence function. It should 
be noted that, in a model in which more than a two-level 
transition is considered, after this time a transition to 
even lower-lying hyperfine spin states are likely. However, 
in our two-level approximation these flips are not taken 
into consideration. 



(42) 



IV. PLANAR MULTILAYER SUBSTRATES 



Up until now, the derivation of all formulas were valid 
for arbitrary substrate geometries. A particular geomet- 
ric arrangement is fixed by defining the correct bound- 
ary conditions for the dyadic Green function G(r, s, w). 
In this section, we will concentrate on the simplest but 
experimentally important realization in terms of planar 
multilayer dielectrics. In what follows, we will focus on 
the spatially nonlocal term in Eq. (|42|) only. In particu- 
lar, we notice that this is equivalent to taking the long- 
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time limit of Eq. I|42|) . Hence, for now we consider only 



S (ri,r 2 ,0M) 



, (Mb9sY 
c 2 eoh 



(i\S q \f)(f\S k \i) 



Im 



V x G{r 2 ,r u u) A ) x V 



kq 



12 



(44) 



which had previously been derived in connection with 
spatial decoherence of matter waves in [3(| . Note that in 
a planar geometry in which the atom is held at a fixed 
distance to the material surface, the spin flip rates Fj 
coincide due to translational invariance, i.e. Ti 2 = T\ = 
IV Note also that Eq. 1(44)1 is temperature- independent. 

Let us first consider a half-space filled with a dielectric 
or metal of dielectric permittivity e(ui) (see the discus- 
sion in Sec. [H}. We evaluate the spin matrix elements 
for the transition from one hyperfine ground state to an- 
other by the basis states through the Clebsch-Gordon 
coefficients \F,m F ) = Y, msmi ^71™^ ™r>- For the 
87 Rb ground state transition |2, 2) — > |2, 1), the non- 
zero matrix elements are |(i|£i/,z|/)| = 1/4. The dyadic 
Green function for such a situation can be found in 
IH H HI- We have collected some of the formu- 
las in Appendix ["""] Note that in the expressions for 
the components of the generalized reflection coefficient, 
Eq. (|A3|) . the common factor e lk iz( z + z ) = e ^k lz d can ^ e 
approximated by e _2d l fe IM because the transition wave- 
length, A = c/(2ttuj), is the by far biggest length scale 
in the system such that the approximation k\ z « — k 2 
holds. Then, by going over to polar co-ordinates in the 
two-dimensional Fourier transform in Eq. (|A2|) . k|| = 
(k x ,k y ) 1— > (K cos ip, K sin ip) and <i 2 k|| 1— > KdKdip. We 
can thus write Ti2 after integration over (p as 



r = {^B9s) 2 

8c 2 £q/i 



K 2 dK e~ 2Kd 

"(27T) 5 2~ 



Im[rf 2 £ ] 



(45) 



It is worth noting at this point that the line shift 8uj a 
in Eq. can be computed as in Eq. (|4~5l) by replacing 
Imfr^] with Ke[rf 2 E ]- Moreover, it is easily seen that 
both T a and 5u> a are of the same order. 

Let us assume that an atom is located at a distance 
d away from the planar interface which we describe by 
its skin depth 5. In our example, we have chosen an alu- 
minium substrate with S — 110/iin and an atomic tran- 
sition frequency as / = 560kHz. Furthermore, the atom 
can be in two distinct positions with a lateral separation 
I. In Fig. n we show the decay of the spatial coherence as 
measured by the function S r 2 , lj a ) for varying sepa- 
ration I in /im for three different atom-surface distances 
d. As a function of separation, the decay of the spa- 
tial coherence starts off rather slowly We attribute this 
behaviour to the fact that for separations below the co- 
herence length of the magnetic-field fluctuations the spin 
flip is driven coherently at both sites. 

In order to investigate the small-separation limit in 
some more detail, we take a closer look at the Weyl 
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FIG. 1: Spatial coherence function of the fluctuating mag- 
netic field S (ri, r2, cja), Eq. (1441 . as a function of the lat- 
eral separation I in /im with the parameters / = 560 kHz, 
S = 110 /im for three different distances from the surface: 
d = 20 /im (solid line), d = 10 /im (dotted line), and d = 5/im 
(dashed line). 



expansion of the scattering Green tensor Ji^ 12 ) (r, r', lu), 
Eq. HA2fl , which is the by far dominant contribution com- 
pared with the free-space Green function. The separation 
I is nothing but I = \g — q'\ and serves as a parame- 
ter in the integral. Hence, we can expand the exponen- 
tial e' k n'( e-e > in Eq. 1|A2() into powers of I and evaluate 
each term seperately. The zeroth-order coefficient triv- 
ially leads to the spin flip rate T^. The contribution from 

terms proportional to I vanish identically due to the sym- 

( 12) 

metry of the generalized reflection coefficients with 
respect to the wave-vector components k|| in the (x, y)- 
plane. In fact, all odd powers of I vanish because of that 
symmetry. 

Hence, the lowest non- vanishing power is I 2 . It is 
straightforward to find analytical expressions for the spa- 
tial coherence in that limit by converting the additional 
factor K 2 coming from the expansion of the exponential 
in Eq. 1A2|) into a parameter differentiation with respect 
to the atom-surface distance d. That is, we make the 
replacement K 2 1— > j -S-^ . In this way we find that 



S (ri,r 2 ,0M) = ( T u 
1 12 



96 dd 2 



0{l A ). (46) 



In certain asymptotic regimes in which Ti2 can be ex- 
pressed as a monomial oc d~ n of the atom-surface dis- 
tance d (see, e.g. 0,0, Eg), Eq. 1 (4*5) 1 can be rewritten 
in the form 



S (t 1 ,t 2 ,uja) = 1 



5n{n+ l)l 2 
96cP 



Oil 4 ). 



(47) 



In addition to the planar half-space we consider the ex- 
perimentally relevant situation in which a thin metallic 
layer of thickness h has been brought onto a dielectric 
substrate. The generalized Fresnel coefficient for this 
three-layer system is given in Eq. I|A5)I . In the limit of 
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thick films (S, h ^>> d) the asymptotic behaviour of the 
spin flip rate is T 12 oc 1/d Il4| whereas for thin films 
(5 > d > h) we have T 12 oc 1/d 2 Hi H^- Thus, we 
finally obtain the small-Z limit of Eq. (|42[) as 

Qi,(t) = l-^(l-e-^)+0(l*), (48) 

where a = 1 for thick films and a = 3 for thin films. It is 
interesting to note that the fall-off is three times faster for 
thin films than for thick films which we attribute to the 
fact that in thick films it is more likely to drive spin-flips 
coherently. 

In order to see how the time scale is related to the ex- 
pected lifetime we can expand the exponential in Eq. Ij48(l 
for short times as 

\ gi2{t )-g 12 (0)\ = ^(^j+O(t 2 ) (49) 

where Qi 2 (0) = 1 and r = T^ 1 . The left-hand side in 
Eq. ([^ can be thought as a proper measure of deco- 
herence due to spin flips in terms of physical parameters 
such as the spin-flip lifetime r, the separation I and the 
distance from the surface d. This means that it is possi- 
ble to maximize those experimental parameters while the 
decoherence rate is under control. Hence, Eq. 1|49[) turns 
out to be particularly interesting from the quantum in- 
formation point of view when a certain degree of spatial 
coherence has to be maintained. 

For larger separations, however, it is difficult to find 
analytical approximations and one has to resort to nu- 
merical evaluations of the Fourier transform l|A2jl . It is 
interesting to see at which separation l\/2, as a function 
of the other length parameters in the system, the spatial 
coherence drops to half its initial value which could be 
taken as a measure of robustness. In Fig. [21 we show the 
dependence of li/ 2 on the thickness h of the intermediate 
layer. In our calculations, we assumed a transition fre- 
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FIG. 2: Lateral separation I1/2 after which spatial coherence 
has dropped to half its initial value as a function of the layer 
thickness h. The skin depth was varied from S — lOO^im 
(solid line) to S = 50fim (dashed line) and 5 = 10/xm (dotted 
line) . The atom-surface distance was d = 50/mi and all other 
parameters as in Fig. Q 



quency of f — 560kHz. We have plotted li/ 2 for three dif- 
ferent skin depths: S = 100/im (solid line, corresponding 
to a good conductor such as Al of Cu at room temper- 
ature), 5 = 50/im (dashed line), and S = 10//m (dotted 
line). Although the latter two skin depth values are not 
realistic for materials at room temperature, at kryogenic 
temperatures these values can be achieved. For example, 
just above its critical temperature of T c — 9.2K, pure nio- 
bium shows a skin depth of only 5 = 15/im at / < 1MHz 

m- 

In Fig. |21 it is clearly seen that for skin depths smaller 
than the atom-surface distance (dotted line), the robust- 
ness of spatial coherence drops dramatically with increas- 
ing substrate thickness h until h ~ S. This can be un- 
derstood when noting that by increasing the thickness of 
the intermediate layer one increases the number of fluc- 
tuating dipoles that can cause the spin flip. Any further 
increase beyond h ~ 5 does not change much because 
fluctuations would not reach the substrate surface. Note 
also that the coherence length l Y j 2 levels out roughly at 
the value of the skin depth, l X j 2 ~ 5. 

For skin depths equal (dashed line in Fig. 0) or larger 
than the atom-surface distance (solid line) spatial coher- 
ence is robust over a wide range of substrate thicknesses 
h. Only for h > 5 does the coherence length decrease 
towards the atom-surface distance. 



V. CONCLUSIONS 

In summary, we have investigated loss of spatial co- 
herence of atomic superpositions due to thermally driven 
spin flips. The consistent quantization of the electro- 
magnetic field in absorbing dielectrics and metals allowed 
us to employ a first-principles approach to decoherence 
in this particularly simple physical system. The quan- 
tization scheme is based on the source-quantity repre- 
sentation of the electromagnetic field in terms of the 
dyadic Green function of the associated classical scat- 
tering problem and a bosonic vector field that serves as 
the dynamical variables of the theory. The Green func- 
tion contains, via the dielectric permittivity, all infor- 
mation about the geometric arrangement and material 
properties of the substrate. Because the theory, starting 
already with Eq. QJ, is strictly valid only for spatially 
locally responding materials, we stress again that spa- 
tially nonlocal effects — which could be non-negligible 
for small skin depths (i.e. large conductivities) and small 
atom-surface distances — have not been considered. 

The interaction dynamics between atomic spin 
and electromagnetic field has been described in the 
Schrodinger picture and the Markov approximation 
which led to the result for the time evolution of the 
off-diagonal matrix element (or coherence) Qi 2 (t) of the 
single-particle density matrix, Eq. 142|) . The spatially 
nonlocal part, Eq. I]44p. agrees with previously obtained 
results [361 ] for spatial decoherence of matter waves. It 
should be noted that both Eqs. and jjljl are valid 
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for arbitrary geometrical arrangements of substrate ma- 
terials. 

For planarly multilayered substrates the dyadic Green 
function is explicitly known |24l IH IHIll and the mam 
formulas presented in Appendix [S] For small lateral sep- 
aration I of the atom's two possible positions we found 
that the spatial coherence decreases quadratically with I 
and inversely proportional to the squared atom-surface 
distance d [Eq. (|48|l ]. For larger separations, a numeri- 
cal study of a three-layer system showed that the coher- 
ence length l\/2, defined to be the separation after which 
the coherence decays to half its initial value, converges 
for thick intermediate layers to roughly the atom-surface 
distance d. 

We believe that these results are important for the de- 
sign of microstructured devices in which spatial coher- 
ences are used to encode quantum information. In par- 
ticular Eq. 149f) shows how the decoherence rate depends 
on experimental parameters such as lifetime, lateral sep- 
aration and atom-surface distance. They can be tuned in 
order to fall within a given tolerance rate for the degree of 
decoherence. Therefore, the theoretical results presented 
here may be useful in the physical realization of atomic 
traps where a certain degree of spatial coherence has to 
be maintained in order to be able to perform some kind 
of error correction. 



parts of the total scattering Green function, respectively. 
These scattering Green functions satisfy the homoge- 
neous Helmholtz equation. In our case, we only need 
to concentrate on the reflection part R^- 12 ^ (r, r', lu). 

The translational invariance in two spatial directions, 
say in the (x, y)-plane, allows one to write the Green 
function in terms of its Weyl expansion 



(27T) 

(A2) 

[p = (a;, y)] where k|| = (k x ,k y ) is the wave- 
vector in the (x, y)-plane. The matrix components of 
I?( 12 ' (k|| , u>; z, z') can be read off from [24| as (here we 
omit the arguments to enhance readability) 



12) 

XX 



2ki 



% _ e ik lz (z+z r ) 



k 2 h 2 k 2 
TM ^lz^x , TE N 

12 T2T2" + r 12 12 



#(12) 



2k! 



1 c ikiAz+z') 



TM 



k{k 
h 2 h h 



12) 



2fci 



1 c ik lz (z+z') 



12 k 2 k 2 

l\<^ hj || 

TM ^lz^x 



k h 

Pl2 ~fcT 



#(12) 
± '•zz 



ikiz(z+z') 



2k 



Iz 



12 k\ 

k 2 

TM II 
12 ~k\ 



(A3) 
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APPENDIX A: GREEN FUNCTION FOR 
PLANAR MULTILAYERS 

We briefly review the calculation of the Green function 
of planar multilayers as it can be found in 24, 37, 38, 39] . 
The dyadic Green function for the electric field scattering 
off a material interface can always be decomposed into 



G(t,v',lu) 



G« (r, r', u) + R {12) (r, r', u) ; r, r' e Vi 
r( 12 )(r,r',w);re Vi.r'e V 2 

(Al) 

where (r, r', to) denotes the solution to the inhomo- 
geneous Helmholtz equation with the source in region Vi 
which in our case is vacuum with £i(w) = 1. The two 
(double-sided transverse) scattering parts i?( 12 ^ (r, r', lu) 
and T^ 12 ) (r, r', to) have to be introduced to satisfy the 
boundary conditions for the electromagnetic fields at the 
interface and describe the reflection and transmission 



where k 2 — ^e^oj) and k 2 z = k 2 — k 2 . The remain- 
ing matrix elements can be deduced from Eq. (|A3|I by 
replacement rules such as Ryy — R^x\k x k y ) and 
the reciprocity condition R^ 12 ^ (r, r', lu) — R^ 12 ^ T (r', r, lu) 
which yields R (12 1 (k,, , lu; z, z') = R^ T (~k l{ , lu; z', z). 

The functions r\® and rf^ 1 denote the usual Fresnel 
reflection coefficients for TE and TM waves, respectively, 
and are defined by 



r 12 



hz - k 2z 
hz + k 2z 



TM 



e 2 (Lu)k lz - £i{Lu)k 



2z 



e 2 {uj)ki z - ei{uj)k 2z 



(A4) 



The Fresnel coefficients obey certain recursion relations 
that permit one to calculate the dyadic Green function 
for arbitrarily multi-layered materials |37l I38L |3jj. In 
particular, the generalized Fresnel coefficient for a three- 
layer geometry reads (for both TE and TM polarizations) 



ri2 



ri2 + r 23 e 



2ik 2z h 



1 — r 2 ir 2 3e 2lk2 * h 



(A5) 



where h is the thickness of the intermediate layer 2. 
This relation has been used in the numerical calculations 
throughout the paper. 
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